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Using the determinant representation for the field-field correlation functions of impenetrable 
anyons at finite temperature obtained in we derive a system of nonlinear partial differential 
equations completely characterizing the correlators. The system is the same as the one for impen- 
etrable bosons but with different initial conditions. The large-distance asymptotic behavior of the 
correlation functions is obtained from the analysis of the Riemann-Hilbert problem associated with 
the system of differential equations. We calculate both the exponential and pre-exponential factors 
in the asymptotics of the field-field correlators. The asymptotics derived in this way agree with 
those of the free fermions and impenetrable bosons in the appropriate limits, k 1 and k — + 0, of 
the statistics parameter k, and coincide with the predictions of the conformal field theory at low 
temperatures. 
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I. INTRODUCTION AND STATEMENT OF RESULTS 

This is the third paper in the series investigating the correlation functions of one-dimensional impenetrable anyons 
at finite temperatures. In the previous two papers [H, Q , we have obtained the anyonic generalization of Lenard's 
formula and determinant representations for the time-, space-, and temperature-dependent correlators of the anyons. 
Here, starting from the obtained determinant representation, we derive the explicit expressions for the large-distance 
asymptotics of the anyonic space correlators at finite temperatures. These expressions demonstrate the crossover 
between the fermionic and bosonic behavior of the correlators when the appropriately-defined statistics parameter 
K g [0,1] changes between the corresponding limits of bosons, k = 0, and fermions, k = 1. The approach we 
use to compute the asymptotics of the correlation functions of impenetrable anyons follows the one developed for 
impenetrable bosons [3, S Q ■ The main technical point of this approach is the large- distance asymptotic analysis 
of the Riemann-Hilbert problem associated with the system of differential equations which can be derived from the 
determinant representation. 

The model of impenetrable one-dimensional (ID) anyons we consider is the limit of infinite interaction strength of 
the Lieb-Liniger gas of anyons with delta-functional repulsion d, 0, H] ■ The fact that impenetrable particles can not 
be directly exchanged in ID systems implies that introduction of the exchange statistics for such particles requires an 
additional convention on the choice of the sign of the statistical phase associated with exchange of each pair of particles 
in the system wavefunctions [T^j . The absence of direct exchanges also implies that local thermodynamic properties 
of N impenetrable ID anyons are independent of statistics in the thermodynamic limit N oo, and coincide with 
those of the free fermions - see, e.g. hj. In the case of quasiperiodic boundary conditions (particles moving on a 
^ [ circle), thermodynamic properties are statistics-dependent at the level of corrections to the leading large- asymptotic 
• ^ . terms El [11 m, [H ■ Although the non-local characteristics like single-particle momentum distribution (i.e., Fourier 
' transform of the single-particle density matrix) depend strongly on the statistics parameter k [l^, [l^ , they may not 
^ , be measurable in practice [l3|- For transparent anyons, i.e. finite particle-particle repulsion, local thermodynamic 
■ - - ' properties do depend on k through modification of the effective coupling constant [ll, Elj El, [III , but the model 
of anyons without a hard core is not well-defined because of the essential singularity of the system wavefunction at 
coincident particle coordinates. Qualitatively, these pro per ties of the one-dimensional anyons are preserved in other 
related anyonic models with, e.g., motion on the lattice ]20l. [2ll. [2^. [23| or different interactions p4l. [25j. 

This work presents the exact calculation of the large- distance asymptotic of the field-field correlator of impenetrable 
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anyons. To state its main result, we need some notations and definitions. First, we introduce the function 

where the brancli of the logarithm is specified by the requirements that no branch cut intersects the real axis, and 



X 



1™ ln(-3^,3^-— ) ^0. (2) 



We then define two constants C(/3, k) and c(/3, k) as 



C(/3, k)^2i I v{\ f3)dX = - / In ( ^^^^ ) dX , (3) 



TT 



and 



c(/3, k) = / — ^ dXdfi . (4) 



OO 



2(A-m) 



In terms of these quantities, our main result for the two leading terms in the large-distance asymptotic of the field-field 
correlator of impenetrable anyons is: 

{^^{xi)^{x2))t = ^C7(VT,«)/2gC(VT,«) |^coe"-^^° + c_ie'^-^^-i) . (5) 

Here T is temperature, h - chemical potential, xi2 = {xi — X2) — > -l-oo, and the constants Xj with j — 0,-1 are 
complex numbers which depend on j3 = h/T and statistics k, as 

A, = (l3+ ^13^ + 7r2[K + 2j]2) /V2 + z (-/3 + V/?' +^^[>^ + '^i?) ^'^ 7^2 . (6) 
Two other constants Cj which give the amplitudes of the asymptotic terms in ([5]) are given by 

_ .e'^"^^/^a-^(Aj) 



A, 



(7) 



where the function a(X) is defined by Eq. (|53p . The second term in the asymptotic expansion ([5|) is always smaller 
than the first one, but becomes relevant when the statistics parameter approaches the fermionic value k = 1. The 
precise conditions on k under which this term does not exceed the accuracy of our calculation are stated in Section 
IV A[ see Eqs. (|72p and ([73]) . More generally, all results of this work are obtained, strictly speaking, only for k e (0, 1]. 
However, as shown in Section IVI A[ in the limit k — > our results reproduce completely the exponential behavior 
of the asymptotics for impenetrable bosons [k = 0) 0, EH, and the pre-exponential factors in the case of negative 
chemical potential h. 

An important role in obtaining the asymptotics ([5]) is played by a set of auxiliary potentials Bj^ and 

introduced below in Section [IT] which are related to the integral operators in the representation of the correlation 
functions. The integral operator which enters in the Fredholm determinant representation for the correlator is of a 
special kind called "integrable integral operator" . Integral operators of this type play an important role in the study 
of correlation functions of integrable models and random matrices fsl. R [29l. [SOl. I3ll.l3^. The special form of the kernels 
of such operators allowed us to obtain a system of partial differential equations for the auxiliary potentials, and for the 
logarithm of the Fredholm determinant, which characterize completely the correlators at any distance. The system 
of equations is the same as the one obtained for impenetrable bosons [1, d, HI] but with different initial conditions. 
The short-distance and low-density behavior of the field- field correlator can be extracted from the initial conditions. 
At zero temperature, the nonlinear differential equation for the logarithm of the determinant becomes the Painleve V 
differential equation obtained by Jimbo, Miwa, Mori and Sato for impenetrable bosons [sot . but again, with different 
initial conditions. For finite system, the Painleve VI differential equation characterizing the field- field correlator was 
derived by Santachiara and Calabrese in [s^ . The short-distance behavior of the 2-point correlation function obtained 
in [s^ l agrees in the thermodynamic limit with our results at zero temperature. The large-distance behavior of the 
correlation functions is extracted from the solution of a matrix Riemann-Hilbert problem associated with the system of 



3 



differential equations for the auxiliary potentials. The technique used is similar to the one developed for impenetrable 
bosons [1, . For the computation of the pre-exponential factors in the asymptotics of the correlators we have 

adopted the method employed by Kitanine, Kozlowski, Maillet. Slavnov and Terras in their study of the generalized 
sine-kernel [36j . This method was first used by Cheianov and Zvonarev in the computation of the asymptotic behavior 
of correlation functions of impenetrable electrons [1^, H^] ■ Some of the results of this work were presented in . 

The plan of the paper is as follows. In Section |lTl we briefly review the determinant representation for correlators 
obtained in our previous papers , and introduce the auxiliary potentials which play a central role in the subsequent 
analysis of the correlation functions. In Section IIIIl we derive the system of partial nonlinear differential equations 
characterizing the auxiliary potentials and the correlators. Scction llVl introduces the matrix Riemann-Hilbert problem 
associated with the obtained system of equations. The large-distance analysis of this problem is performed in Section 
IVl We conclude in Section IVII by comparing our main result with the known limiting cases: the low-temperature 
limit characterized by the conformal behavior of the correlators, and bosonic and fermionic limits of the statistics 
parameter. Technical details of the calculations are presented in the appendices. 

II. DETERMINANT REPRESENTATION AND AUXILIARY POTENTIALS 

The model of one-dimensional anyons we consider in this work is characterized by the hamiltonian H given in the 
second quantized form as 

H= ( dx {[d.^'^\x)][d^'^[x)]+c'^\x)'^\x)^[x)^[x) - h'^\x)'^{x)) . (8) 



Here c is the coupling constant which is assumed to be large, c —t +oo, so that the anyons are impenetrable. The 
fields in the Hamiltonian ^ obey the anyonic commutation relations 

^'(a;i)*^(a;2) = e-''"'<'-'^-^^''-^\x2)-^{xi) + 5{xi - xa) , 

where eix) = xl\x\, e(0) = 0. The commutation relations become bosonic for k = 0, and fermionic for k = 1. In this 
work, we are interested in the large-distance behavior of the field-field correlator of the impenetrable anyons ^ . The 
static, i.e. equal-time, correlator at finite temperatures is defined by the standard relation: 



(*T(xi)*(x2))t = Tre-g/^ ■ 

In the previous paper of the scries, we have obtained the following expressions for this correlator 



(9) 



{^^\Xi)^{X2))t = 

and 



*^(a;i)«'(x2))T = ^Tr (1 - 7/^t)"'A+ det(l - 7AT)|^=(i+e+"»)/. , xi > X2 , (10) 



{¥{x^)^{x2))t - ^Tr [(1 - ^KtY^A-^ dct(l - 7/^T)^ = (l+e— )/. ' < X2 . (11) 

In these expressions, Kt and are the integral operators acting on the entire real axis with the kernels 

A — jJL 

where 

^(A)^^(A,T,/.) = ^^-^l^, (12) 

is the Fermi distribution function, and Tr[/(a;,y)] = J f{x,x) dx. Equations pop and pip remain valid at zero 
temperature, but in this case, the integral operators act on the interval [— q, q] with q = \/h, and have kernels 

if (A,;,) = , A±(A,/i) = eT"-(^+'^) . 

A — /i 
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Equations (fTO|) and ((TT|) show that the anyonic corrrelator (^''''(xi)4'(x2))t depends on the sign of ii — 12. However, 
since the values (jlOp and (jlip for positive and negative Xi — X2 are related directly via complex conjugation, one can 
focus only on one range, e.g., the correlator pO|) . This correlator depends on four variables: the coordinate difference 
xi — X2 > 0, temperature T, chemical potential h, and statistics parameter n. As we will see bellow, introducing the 
variables that are rescaled by temperature: the distance x and the chemical potential (3 defined as 

X = ^{Xl - X2)VT , P=j,, 

and similarly changing the spectral parameter, A XVT, one makes explicit dependence on temperature very simple: 

(vl/t(.Ti)*(x2))T = ^ff(a:,/?,7)l7 = (l+e-)A- (13) 

The function g{x,P,"f) here will be defined in the next section. 



A. Auxiliary Potentials 

The Fredholm integral operator Kt appearing in the expressions (jTU]) and (fTTI) for the field correlators belongs to 
a special class of "integrable" operators [1, 0, [121 ■ This means that in terms of the "plane waves" introduced as 



e±(A) = v^e±^^^ (14) 
the kernel of Kt can be written in the following form which generalizes the simple factorizable kernels: 

/X N e+(A)e-(M) - e-(A)e+(/x) 

^^(^''^^^ • ^ ^ 

This is a particular case of the more general situation studied in 0, [lil (sec also Chap XIV of 0]). An important 
feature of this class of integrable operators is the fact that the kernel -Rt(A, ^j.) of the resolvent operator Rt, defined 
by the relation 

Rt^{1~jKt)~^Kt, i.e., (l-7/<T)(l+7i?T) = 1, (16) 
has the same form 1,13. Indeed, Eq. dUl) means that the resolvent kernel Rt{X,^) solves the integral equation 

/•+CO 

RT{\lJi)-l / KT{\v)RT{y.^l) dv = KT{\^i). 



Introducing then the functions /±(A) which are the solutions of the similar integral equations 

p + CO 

f±W-il KT{\^l)f±{^l)d^l^e±{\), (i7) 

T'- 



one can show (for a proof, see Chap XIV of Q ) that the resolvent kernel can be written in the same form (fT5|) as K' 



u ^ /+(A)/-(aO-./-(A)/+(a^) 

An important role in our asymptotic analysis of the correlator pop is played by the auxiliary potentials Bim which 
are defined by 

r+oo 

Bi,„{x,(3,K) = -f ei(A)/™(A) dA, I = ±, m = ± , (19) 



where 7 = (1 + e^'^'^)/Ti. A new feature of the auxiliary potentials Bim in the case of anyons in comparison to bosons, 
when 7 = 2/7r, is that they are now complex. The potentials, however, still satisfy the relations {xJ3,k) = 
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B ^^{x, P, k) and B+j^{x, (3, k) = B (x, /?, k) as in the bosonic case. Indeed, as one can see from the definitions (fT9|) 

and (fT7)) 

B+- = 1 e+(A)/_(A) dA = 7 / e+(A) / (1 ~ 7At)-'(A, M)e-(M) 



~00 J —OQ 

f +00 

7/ e-(A)/+(A) dA = i?-+, 



where in the last Une we have used the fact that the kernel Kt is symmetric, Kt{\, fJ-) — KxifJ-, A). In order to prove 
the second assertion we start with the integral equation p7|) for /+(— A) 



+ 00 



/+(-A) - 7 / KTi-x,^l)f+{^i)d^, = ^JiXje' 

which can be rewritten as 



+00 

iXx 



/+(-A) - 7 / A't(-A, -fi)f+{-fi)dfi = ^Mxje 
<y ~Qc 

Using the relation Kt{~X,^IJ.) = Kt{X,ij.), we see then that /+(— A) = /-(A). This gives 

/+00 r+oo 
e+(A)/+(A) dA = 7 / e+(-A)/+(-A) dX 

/-f oo 
e_(A)/_(A) dX = . 
-oo 

Finally, the definition of the kernel shows directly that pO|) can be rewritten as 

(vI't(xi)^(x2))T = ^det(I-7/^T)U=(i+e+"»)A / /-(A)e-(A) dA, 
in terms of the function /-(A) pT|) . This shows that the function g{x,f3,j) introduced in Eq. ((T3)) is given by 

ff(a^,/3,7) = B++{x,(3,j)dct{l ~ jkT)\j=(i+e+i'''}/7r ■ (20) 

III. DIFFERENTIAL EQUATIONS FOR STATIC CORRELATORS 

In general, it is difhcult to obtain directly differential equations for the entire correlator (|10p . The strategy we 

pursue is to obtain first nonlinear partial differential equations for the potentials , B^ , and then show that the 

function (t{x,(3,j) = lndet(l — jKt) can be expressed in terms of B++ and _B_| We start by look for the two 

operators L and M depending on -B++ and B^ such that we have 

d,F{X) = L^^(A), 
{2Xdp + dx)F{X) = M ^^(A) , (21) 



where F(X) is the two-component vector function 



/+(A) 



^^(A) = ^ ) . (22) 

The compatibility condition for the system (PT|) 

(2Aa^ + 9A)L~a,M + [L,M] =0, 

gives then a system of partial differential equations in x and /3. It should be noted that since the Fredholm determinant 
in Eq. (jlOp differs from the one appearing in the similar representation for impenetrable bosons only in the value of 
7, the calculations follow closely those for impenetrable bosons [1, [s^. The main difference is that 7 for anyons is 
complex, making the auxiliary potentials -B++, -B-i complex and not real as for the bosons. 
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A. The L operator 



We start from the integral equations for the functions /±(A) 

f +00 



/±(A)-7/ KT{X,^i)f±i^i)d^l = e±iX), j = {l + e^^-)/n. 
Differentiation with respect to x gives 

(f - 7^T)a,/±] (A) - 7 [dJ<Tf±] (A) = ±i\e±{\) . 

Acting on both sides with (1— 7AV)~^ and using the relations (1— 7AV)^"'^ = (1+7-Rt), (1 — jKT)^^e± (A) = /±(A), 
dxKriX, fi) = (e+(A)e_(/^) + e_(A)e+(/x))/2, and the special factorization for the resolvent kernel ^TE\i wc obtain 

a,/+(A) = *A/+(A)+i3++/_(A), 
5,/_(A) = -zA/+(A) + i3__/_(A). 

These equations mean that the L operator has the form 

L = ^Aa3 + B++ai , (23) 

where we have used the fact that = B , and ai are the Pauli matrices 



CTl = 



1 

1 



0-2 = 



-i 

1 



CT3 = 



1 

-1 



B. The M operator 



The derivation of the M operator is more complicated. In this case, we have to rely heavily on the following property 
of the Fermi distribution function ^?(A): 

{2Xdp + dx)^{X) ^ . (24) 

This property is essential, since in the derivation of the M operator, one encounters terms which contain the resolvent 
Rt{X,ij,) but still can not be reduced to the form (A — fi)RT{X,fi) (as in the previous section) which produces 
"disentangled" terms (i.e., products of the one-dimensional projectors). This is also why the differential operator 
associated with M is (2A9/3 + d\) and not simply djj as one could have expected. 

The computations are more involved but similar with the ones for impenetrable bosons and we refer the reader to 
Chap. XIV of [aj for the derivation. The final result is 

{2Xdi3+dx)f+{X) = txf+{X)-if+{X)dpB^++tf^{X)di3B++, 
{2Xd^ + 9a)./- (A) = -ixf+{X) - if+{X)dpB^^ + tf^{X)dpB+_ . 

Taking into account that B |- = , -B++ = B this finally gives the M operator as 

M = ixa^ — idpB^ 0-3 — dpB++(T2 ■ (25) 



C. Differential equations for the potentials 

The results obtained in the previous sections allow us to state the following theorem: 

Theorem III. 1. Forallj — (l+e*'^'')/7r with k ^ [0:1); the potentials B^ {x, (3,j) , B++{x, (3,j) satisfy the following 

system of partial differential equations 

dpB+^ = x+- — , (26) 

2 B++ 

dxB+^ = Bl+, (27) 
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with the initial conditions ( at fixed (3 ) 

B++(x,/3,7) = i?+-(a;,/3,7)=7rf(/3) + M(/3)]'a^ + 0(a;2), x^O, (28) 

where d{(3) = •d{X)dX and 

B++ix,p,j)^B+^{x,(3,j) = 0, (29) 
The potential i?++(x,/3, 7) satisfies the nonlinear equation 

with the same initial conditions. 

Proof. In the previous sections we have shown that the two-component vector function F{X) satisfies the following 
differential equations 

9,F(A) = LF(A), 
{2X8(3 + dx)F{X) = M Fix). 

where L,M are given by and The compatibility condition for these equation is 

[a, - L, (2X8(3 + dx) - M] = {2Xdf3 + dx)l - d.,M + [L, M] = . (31) 

Differentiating Eqs. (l23|) and ^E^, we get 

i2Xdp + dx)L = 2XidpB++)ai + taa , 

dxM = i(T3 - i{dxdpB+^)a3 - {dxdi3B++)a2 . 
Using the standard relation for the Pauli matrices 

[cTi, CTj] = 2i£ijk(yk , 

we also have 

[L, M] = -2X{dpB++)ai + [2xB++ - 2B++{dpB+^)] - 2iB++{df3B++)cj^ . 



Plugging these results into Eq. ((3T|) we obtain a matrix with only two distinct elements. The condition that this 
matrix vanishes identically gives the following relations 

2B++{di3B+^) - 2xB++ + d^dpB++ = , 

and 

d^dpB+^ - dpBl+ = . 

The first relation is just Eq. while the second is equivalent to Eq. (P7|) due to The equation ([50]) is obtained 
by differentiating (|26p with respect to x, (|?7|) with respect to /?, and then equating the right-hand sides of the resulting 
relations. The asymptotic behavior at short distances (|28p and at low density (p9)) . which give the initial conditions, 
are obtained in Appendices and (|B]). □ 

D. Differential equations for a 

Theorem III. 2. For any 7 = (1 -f- e'"^'')/7r with k G [0, 1), the partial derivatives of (j{x, f3, 7) = lndet(l — jKt) with 
respect to x and (5 are given by 

d^a = , dla = -B\+ , (32) 

dpa = -xdpB+^ + \{dpB+^f - \{dpB++f . (33) 
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Furthermore, for all js, the function cr(x, /?, 7) satisfies the following nonlinear partial differential equation 

{dfidlaf + A{dla)[2xdf3d,a + {dpd^af - 2dfia] = 

with the initial conditions 



(34) 



-^d{P)x - [-fd{l3)f 



-00 , 



(35) 



where, as before, d{f3) = 'd{X)d\. 

Proof. We start with the simpler case of derivative with respect to x. Making use of the following standard represen- 
tation for the Fredholm determinant of an arbitrary integral operator 



det(l - jKt) = exp I - V — Tr if? 

1 ^-^ n 



n=l 



(36) 



where Tr Kt = / Kt{\, A), Tr = J Kt{X, ^)Kt{ij., A) dXdji, and so on, we obtain 



d.,a = -7 Tr (1 - -ikT)'^d.,kT 



From the relations 
used before, we have 



(1 — "fkr) ^e± (A) = /±(A) and dxKT{X,l^) — (e+(A)e_ (/it) + e_(A)e+(/x))/2 that were already 



-7 Tr 



(1 - jKry^d^kr 



+ 00 



/+(A)e_(A)dA 



+ 00 



/_(A)e+(A)dA, 



= --(i3„+ + i?+_). 



which due to B |- = means that 



proving the first part of Eq. (|32|) . The second part is obtained differentiating once more with respect to x and using 
Eq. 123). 

The derivation of the differential equation ()33p involving the P derivative is more cumbersome. Again, the compu- 
tations are similar with the impenetrable bosons case and we r efer the reader to Chap. XIV of 3. 

The last differential equation ([34]) is obtained by replacing the potentials with B^ = —dxO' and B++ = (— O^cr)^/^ 

in the R.H.S. of Eq. ([55]) . It is also straightforward to see that KT{X,fi) — > when x ^ or /3 — > ~oo, which means 
that (7 = lndet(l — jKt) = in these limits. The first condition in ([55]) follows then from Eq. (pS)) . This concludes 
the proof. □ 

Finally, we end this section with a simple but important observation. Integrating equation dxC = —B^ with initial 

conditions l[35|). we express the function g(x,/3,7) (|20|) and, therefore, the field correlator ^3]). directly in terms of 
the potentials , B^ : 



g{x, /3, 7) = B++ [x, /3, 7)e- -^0^ ^■'5,^)^, 
This expression will be used in our subsequent analysis of the field correlator. 



(37) 



E. The zero-temperature limit 

In the zero temperature limit, the field-field correlator (|13p depends only on three variables, the distance (xi —X2) > 
0, the Fermi momentum q = \fh (or, equivalently, the chemical potential /i), and the statistics parameter k. The 
distance and momentum dependence can be encoded by one variable: 

C = ^^l^q = xP^I- , 
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where now the rescaled variables are x = (a;i — X2)/2 and f3 = h. At T = 0, the integral operator K acts on the 
interval [—q, q] and has the kernel K{X, n) = sin a;(A — /i)/(A — /i) . The logarithm of the determinant is 

o'o(^,7) = lndet(l - 77f)|^^(i+e""')/7r • 

The partial differential equation (j34p characterizing a(x,f3,j) becomes then an ordinary differential equation in ^. 
Introducing the new function 

<Jo = a^oY (38) 
where prime denotes the derivative with respect to ^. we transform Eq. (j34p into 

{^<yo? + 4(e'^o - cToMao - 4^0 + {a'„f] = (39) 

with the boundary conditions 

ao = -27^-472^2^0(^3)^ 7=(l + e*--)/7r. (40) 

The ordinary differential equation p9p is the same Painlevc V equation obtained by Jimbo, Miwa, Mori and Sato in 
their celebrated work on the one-particle reduced density matrix (field- field correlator) of impenetrable bosons [30| . 
The only difference we find in the anyonic case is the boundary conditions (|40p which, unlike the differential equation 
itself, depend on the statistics parameter. In a certain sense, this result could be expected, since it was noticed 
already by Jimbo et al. that the same equation (j39p . but with different boundary conditions, characterizes the density 
matrix of both the imp enetrable bosons and the free fermions. Similar situation was also noted by Forrester, Frankel, 
Garoni and Witte [35j in their study of systems of finite number of particles with periodic boundary conditions. In 
this case, the reduced density matrix satisfies a Painleve VI differential equation. Their work was recently extended 
to impenetrable anyons by Santachiara and Calabrese [s^l- At T = 0, our results agree with those of Santachiara 
and Calabrese. Indeed, Eq. (|A5p for the first terms of the short-distance expansion of the field-field correlator, which 
reflects the boundary conditions (j40|) . reduces at T = to: 



{¥{xi)^ix2))o = Do(^l- ^Dl{x, - x^f + l^Dl{x^ - X2f^ + O ((xi - x^f) , 7 = (1 + e^^h , (41) 

where Dp = q/ir is the density at T = 0. This result coincides with the expansion obtained in [s^] - see Eq. (47) of 
that work and note that the definition of the statistical parameter k there differs in sign from ours. 



IV. THE MATRIX RIEMANN-HILBERT PROBLEM FOR THE FIELD-FIELD CORELATOR 

The most difficult step in the analysis of the field-field correlation function is the calculation of its large-distance 
asymptotics. As we have seen in the previous sections, the correlation function is characterized by a completely 
integrable system of nonlinear partial differential equations - see Thm. Illl.li (|13p and ((37|). A powerful method of 
analyzing these differential equations is the matrix Riemann-Hilbert problem (RHP) formalism [38|. The solution of 

the associated matrix RHP will allow us to obtain the large-distance asymptotics of the potentials ,^+-(-, and 

therefore, the large-distance asymptotics of the correlator. 

We start by formulating the matrix Riemann-Hilbert problem relevant to the case of potentials , In what 

follows, K G (0, 1]. As we will see below, in this case we need to find a 2 x 2 matrix- valued function x(A), nonsingular 
for all A G C and analytic separately in the upper and lower half-plane, which is equal to the unit matrix at A = 00 

x(oo) = / , 

and has the boundary values on the real axis which satisfy the condition 

X_(A) ^ x+(A)G(A), x±(A)= lim Y(A±*e) AeR. (42) 

The matrix G(A) is called the conjugation matrix associated with the RHP and is defined only for real A. In our case, 
it has the form 

G(A) = f 1 + "^^+i^]f- ' , ""^f (43) 
\ 7r7ei(A), 1 - 7r7e+(A)e_(A) J ^ ' 

where all notations are the same as before, e.g., e-i-(A) = '&{X)e^'^^^ , etc. This means that both G(A) and x(A) 
depend also on x, /9 and k, but this dependence will be suppressed in our notations. 
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A. Relation to the auxiliary potentials 

111 this section, we show how the auxihary potentials S-i , S++ are obtained from the sohition of the RHP (|42p . 

The normahzation condition x(oo) = / means that the solution of the RHP has the following expansion for large A 

x{X)^I+^ + o(^), (44) 



A \X\ 

where ^i{x,f3) is a 2 x 2 matrix which depends only on /3 and x. Then, the following symmetry of the conjugation 
matrix 

G(A) =aiG-i(-A)ai, (45) 
that can be checked explicitly from the definition (|43p . determines the structure of ^'i. This matrix can be written as 

where the factor l/2i is introduced for convenience, and at the moment we do not know yet that and in 

Eq. are the auxiliary potentials p^ . To prove this, we use the following formulation of the RHP. One can show 
directly (see, e.g. Chap. XV of 0]) that the matrix RHP is equivalent to the system of singular integral equations: 

lp^±Mi_^.,, ACM. (47, 

27ri ^ - A - lO 

If we define 

x(A)=x+(A)i?(A), (48) 

where E{X) is the triangular matrix 

1 e+(A) 
e_(A) 



the system of integral equations (j47p can be rewritten as 



with 



X(A) = EiX) + -L r m^d, , A e M , (49) 
Zm 7_oo /-t — A — «0 

GiX,,)=E-\,)[I-Gi,)]EiX)=^_^^l^^^ .,(e,(;.)e_(A)'-e,(A)e_(A^)) 
For the matrix elements, Eq. gives the integral equations: 

/+00 
-oo 

/ + 00 
-oo 

with the same kernel we considered before: 

e+(A)e-(pt) - e-(A)e+(/i) 

= ' 

and 

Xii A) = 1-777/ r TT^/^' 

X2l(A) = -— / r^dfi. 
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In the first two equations, the nonsingular character of the kernel Kt makes it possible to neglect the infinitesimal 
shift in the denominator. These equations are just the integral equations (|17p defining /±(A). Therefore, we can make 
the following identification 

Xi2(A) = /+(A), X22(A) = /-(A). (50) 

Also, noting from Eq. (|48p that xii('^) = Xii,+ (-^) ^-J^id X2i(A) = X22.+ (A), and are analytical in the upper half plane 
of A, we have 

Xii(A = l-^TT / r d/i, 5A>0, 

Taking the limit A — > oo in these equations, we obtain 



Xii(A) = 1 + ^|_ "/+(M)e-(A*)d/i + 0(^^), 
X2i(A) = J^y"'"/_(^)e,_(^)d/i + 0(^^^ 



which shows that the components of the matrix ^'i are indeed the potentials p9|) . as assumed in Eq. (|46|) . As a 
byproduct of this derivation, one can also see from Eq. ((50|) that a part of the matrix product ([48]) means that 

/_(A)j-^+(^)(,e_(A) 
This relation will be used in Sec.|V]in the analysis of the large-distance asymptotics. 



B. An useful transformation of the RHP 



As we have shown above, the potentials B+^,B^ that completely characterize the field-field correlator ([13]) 

through Eq. (|37|) . can be obtained from the expansion of the solution of the matrix RHP ([1^ . This means that the 
large-distance asymptotics of the potentials, and hence of the field correlator, can be obtained from the solution of the 
RHP. To find this solution, we will use equivalent formulation of the RHP in terms of the singular integral equation 
(j47p . As a first step, we transform the matrix RHP (|42p into a new form which has conjugation matrix with 1 on the 
diagonal. This is achieved by introducing the matrix function <i> defined as 

.(A,^.(A,('-r)J,„), 

where the functions a(A),/?(A) are: 
and 

The branch of the logarithm in these equations is specified by the requirement that In (...) for /i — > cx), and it is 
important at this point that k e (0, 1]. Both functions a(A) and /3(A) are analytic separately in the upper and lower 
half-plane, and are the solutions of the following scalar Riemann-Hilbert problems (see Appendix \D\i 

a_(A) =a+(A)g<,(A), Ae R; a(oo) = 1 ; ga(A) ^ + 1 ' ^^"^^ 
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and 

^_(A)=/3+(A)g;3(A), AeR; /3(oo) = 1 ; g^A) ^ e^-^-/^ + 1 ' ^^^^ 
They also have the property 

a-i(A) a(-A), /?-^(A) = /?(-A) . 
Using these relations one can see from Eq. (|42p that the matrix <I> satisfies the following RHP 

$_(A) = $+(A)G<[.(A) , AeM; $(oo) = 1 , (56) 

with the conjugation matrix 

/ 1 -^7^(A)/?+(A)a=i(A)e2^^- \ 

" 1^^779(A)a+(A)/3zi(A)e-2»^^^ 1 ^^^^ 

Equations (|44p and P5|) . together with the fact that q;(A) , /3(A) ^ 1 at A ^ oo, show that in terms of the components 
of the matrix <f> that solves the transformed RHP the potentials B ^ are expressed as 

B+_ = 2i lim A [$ii(A) - P^\X)] , B++ ^ ~2i lim A$i2(A) . (58) 

A— *C50 A— >oo 

The RHP ([5^ and formulae ([55)) will be the basis for our analysis of the large-distance asymptotics of the field-field 
correlator of impenetrable anyons. 

V. LARGE-DISTANCE ASYMPTOTIC ANALYSIS 

In this section, we perform the large-distance asymptotic analysis of the RHP ([551) • The strategy is to use Eq. ([55]) 

to obtain the auxiliary potentials B^-^.{x, /?, k) and _B_| (x, /?, k) from the large-A expansion of $11^+ (A) and $12,+ (A). 

Then, the function a{x,P,K) can be calculated using the differential equations ([5^ and ([55)) . 

A. The auxiliary potentials in the large- a:: limit 



The fact that the conjugation matrix G$ of the RHP ()56p has 1 on the diagonal simplifies the matrix structure of 
the kernel / — of the integral equation formulation ([TTl) of this problem. Combined with the explicit form ([57)) 
of the off-diagonal elements of G$, and Eqs. ([M)) and ([55)) for a and /3, the matrix integral equation gives then the 
following equations for $ii^+(A), $i2,+ (A): 

* fy^ 1 (1 + '^"') '^12,+ (m) «-(m) e-^'^- ^ 

*"'+^^^ = '"^^y-oo (a.-A-zO)/3-(m) (eM--/^-e-) ' = ° ' ^^'^ 

and 

*^^'+^^^ = ^^y_ (M-A-zO)a4M) (eM-^/^-e-) ' = ° ' ^''^ 

An important feature of these equations is that the analytical properties of a(A) and /3(A) together with the fact 
that $11(00) = a(oo) ~ /3(oo) = 1. make it possible to obtain an estimate of $12,+ by shifting the contour in the 
upper half-plane and evaluating the integral by the sum of the residues. Strictly speaking, this approach requires the 
functions <f> to conform to estimates of the type |3>i2(A, x,P)\ < G/A and |$ii(A, /3, a;) — 1| < D/X for QA > in some 
ranges of /3 and x: P < Pq , x > xq, where G, D depends only on the constants Pq ,xo- In what follows, we assume that 
these estimates hold, noting that they can be justified self-consistently as in Sec. 4 of The poles of the integrand 
in Eq. ([50)) are given by A + iO and the zeros A^ of the function e^~~^ — e*'^" in the upper half-plane. The zeros of 
this function are given by the formulae: 

(5iAfe)2 = ^ (/3 + x//32 + ^2[^ + 2fc]2) , (9Afc)' = ^ + ^/WT^^^W+W) , k^O, ±1, ±2, • • • , (61) 
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and, as we will see below, the leading terms in the asymptotic behavior are determined by and A^j^ which are 
closest to the real axis: 

A+ = (/3 + v//32 + 712^2^ / V2 + z (-/3 + + TT^K^) /\/2 , (62) 

and 

All V/3'+^'['t-2]2) /\/2 + i (-/? + V/32 + 7r2[At-2]2) /V2 . (63) 

The superscript + denotes the zeros in the upper half-plane. Closing the contour in this half-plane, we obtain 

$i,,+ (A) (1 + e-)a>n,4A) f;^|^|(^,.j;"'V) + ^^^^^ ' ("'^ 

with 

Q+QA^ (l + e^"-)/3(A+) j>n(A+)e2'^t-- 

^ ^ 2e^-'^ aiK) (A^-A)A+ • 

The series S'^{X) is uniformly convergent for A e M, oiq < < k < 1, and with (3 < Pq. Substituting Eq. into 
([59|) . and using the relations (|54p and ([55|) . we obtain the following representation for Eq. (|59p : 



27ri //-A-zO(eA''-/3-e»'^«)(eA'"~/3 + e«'^« + 2) 

1 /■+°° i?+(/^) 



2ni J_ao A* — A — iO 



d/^ , (65) 



where 



i?+(/i) = (1 + e"") ^"/^ ^ , ^5+(m) . 

The form of the integral equation (p5|) suggests (see Eq. (jPSP in Appendix |D]) that it is useful to consider the following 
inhomogeneous scalar RH problem for the function $(A): 

$_(A) = l>+(A).g(A) +i?+(A), AeK, $(oo) = 1 , (66) 

such that $+(A) = $11. + (A) for A G M. The function g in this RHP is given by the relation 

. (1 + e— )2 

which shows that using Eqs. and ([55)l . p(A) can be written as 

' (e^'-'3-e''^'*)(e^'-/3 + e*^« + 2) a;i(A)/3+^(A) ' 
This form of g{X) combined with Eq. (jD4p implies that the solution of the RH problem ((66l) is: 

^(A) ^ a-^(A)r^(A) - ^-'^'^^-'^'^ r ^-M/^-M^M,,^ ,,c/R. (67) 

27ri ^ - A 

The functions $(A) and $11 (A) are analytic in the upper half-plane. They have the same behavior at infinity: 
$(00) = $11(00) = 1, and their boundary values at the real axis $+(A), 'i>ii_+(A) , A S M are equal. This means that 
they coincide, $11 (A) = $(A), for 9A > 0. Therefore, from Eq. (p7)) and the explicit expression for -R^(A) we have 



$n(A) = a-i(A)/3-i(A)-^^^^^ 2 ^(A^A)$n(A+), 3A>0, 



k— — oo 
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where 

For a; ^ oo, the integral in the last equation can be estimated as C\e~'^^^o where Aq is the k = zero (|6ip in the 
lower half-plane, so the leading term of $11 in the limit a; — > 00 is given by 

4'ii(A) = a-i(A)/3-i(A) + 0(e~'''-*^(^^), 3A > , (68) 

where we took into account that the magnitudes of the imaginary parts of A;^ are the same. Using Eq. (|S5|) in (|60p . 
we obtain the leading term of $12 for large x: 

(;,-A-zO) (eM-^-e-)+^(^ ) , > . (69) 

Substitution of the asymptotics (p5|) and ([55)1 into Eq. ([55)1 gives the large-distance asymptotic behavior of the 
potentials B^+.B^ : 

= - In rf/i + O(e--o ^) , , ^ 00 , (70) 

where the branch of the logarithm is fixed by the requirement that ln(...) for /i ^ 00, and 

The expression for -B++ can be made more precise by closing the contour in the upper half-plane and evaluating the 
integral as the sum of the residues 

max _2 / \ -t- \ 

B++ ^ i{l + e-'^'') ^J^e^'^^ + 0(6-4^^0") , (71) 

k—niin ^ 

where the accuracy of our calculation of the leading term in the asymptotics implies that the sum over poles can be 
limited to the interval between k = m.in[(3, k) < and k = max{P, k) > that are defined by the relations 

^^min < 23Ag , 'iX^n^n^i > , (72) 

and 

SA+,, <29A+, 3A+_+i>29A+. (73) 

B. Determination of a{x,P,K,) in the large-a- limit 

The large-distance asymptotics of the auxiliary potentials obtained in the previous section and the differential 
equations ((32)) and ((33|) give the following asymptotic behavior of the hmction (j{x,/3,j) = lndet(l — jKt)- 

a{x, f3, k) = -xC{(3, k) + c{f3, k) + 0(e-'*'-*^^ ^) , x ^ oo , (74) 

where C(/3, k) is defined by Eq. and c(/3, k) is a constant that depends on /3 and k, but is independent of x. In 
order to determine this constant, which through cr(x, /3, k) determines the (3- and K-dependcncc of the amplitudes of 
the asymptotic terms in the expansion of the field correlators, we integrate the relation 

a^cr = - / i?T(A, A) dX, (75) 

which is obtained by taking the derivative of Eq. ([351) with respect to 7 and using the relation (1 — '^Kt)^^Kt = Rt- 
As shown in Sec. UTTCl Rt{X, A) in Eq. (|75l) is: 

i?T(A, A) = i(9A/+(A)/_(A) - /+(A)aA/-(A)) . (76) 
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Our approach is similar to the one employed in Ref. [3^1 (see also [201), where the authors considered a generalized 
sine-kernel at zero temperature and studied the large- distance asymptotic behavior of the Fredholm determinant using 
the nonlinear steepest descent method of Deift and Zhou [13]. At zero temperature, dct(l — 7^) becomes a particular 
case of the Fredholm determinant considered in [36| . 

The first step is the computation of /±(A) using the relation l(52|) and the following expression for the function x{X) 
in this relation: 



X(A) ^ <i>(A) g 



/3(A) 
«(A) 



(77) 



In the previous Section, we have obtained Eqs. ([55)1 and ([55)1 for the matrix elements $11 (A) and $12 (A) of the matrix 
'I'(A) in the large- .t limit. Closing the integration contour in the upper half plane in Eq. (|69p . one gets more explicitly: 



One needs to find similar asymptotics for $21 and $22- As in the previous Section, the integral formulation (|47p of 
the RHP gives the following equations for these matrix elements: 



1 [ °" *22.+ (m) "-(m) 



and 

$22,+ (A)-l+2-y_^ ^_A-zOa+(M) (eM-^-e-)^''' ^^^^ 
In the integral equation (j79p . one can close the contour in the lower half plane obtaining 

<f>22,+ (A^) a^{\-)e-^\^^-^\^ 



$2i.+ (A) = -7r7 ^ 



k=-c 



A--A-zO/3_(A-) 2A^e— ' 



where A^. are the zeroes of e^" ^ — e"** in the lower half plane. Substituting this result in the integral equation for 
$22,+ and closing the contour in the upper half plane we find 

$22,-h(A) = 1 + 0(e-2(|aA-|+aA+)x) ^ (g^^ 

and, combining Eqs. (jSTjl and ((79)l . 

$2i.+ (A) =0(e-2|'-*VI-). (82) 

Now we are able to compute /±(A) up to exponentially small corrections in x. Explicitly, Eqs. ([5^ and ([77]) give: 

/+(A) \ _ /$n,+ (A)/3+(A)e+(A) + <i>i2,+ (A)a+(A)e_(A) 
/_(A) j - U2,_+(A)/3+(A)e+(A) + <i>22,+ (A)a+(A)e_(A) 



and, therefore, using Eqs. dSH]), dTS]), (gl]), and (EH) we find: 

(83) 



/+(A)\ _ /a;i(A)e+(A)e2-MA) 
/-(A) j ~ "+(A)e-(A) 

The function v(\) here (we have suppressed the dependence on /3) is defined in Eq. ([1]). One can see directly that 
v{X) can be written in terms of the rescaled Fermi weig ht e{\) = (1 + e^ -^)-i as 

KA) = -:^log(l-^7^(A)) = — log' ^ ^ 



and also that 



a+(A) = exp { iTTiy{X) + P.V. / ^^d^ 



+ 00 

/i — A 
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Using Eqs. ^ and ^ in Eq. dTS]), we have 



-\-oo p~\-00 



— OO 

d^u{\){2ix-2dx\oga+{\) + 2iTTdxv{\)) dX. (84) 

— OO 



Our current task is to write the RHS of Eq. (|84p as a derivative with respect to 7. The term depending on x already 
has this form, but the last two, a;-independent, terms do not. We transform them as follows: 



— 00 

-\-oo 



/i — A — iO /i — A + iO 



—00 



^ ' (A - ^ - z0)2 (A - /i + i0)2 J 



Next, we show that this expression is equal to 

2{X-fi) 



I dxyjXyiii) - v{X)d,{y^) ^^^^ ^ 



Indeed, 



1 r+°° 1 1 



— OQ 



(X-n-iO)^ (X-fi + iOfJ 



Therefore. 



d^a = ~2ixd-i / u{X)dX + ^ ' !; ^ dXd^i . 



2(A-m) 

Integrating this relation with respect to 7 and taking into account that (7(7 = 0) = 0, we obtain 

a[x, p, = ~xC{p, K) / — dXdfi. (85) 



C50 



2(A-m) 



Comparison with Eq. ([7^ shows finally that 



c(/3, '^j = y - fj,) ^ ■ ^ ^ 

One can also obtain an alternative expression for the constant c(/3, k) by a simpler calculation. Using again the 
differential equation ((33l) and expressions for the potentials (|70p and ((7T|) . we have 

/?, k) = ^) + lf (^^^^7^) ^Z?' + c(^) + 0(e-4'-^o+-) , ^ ^ oo , 

where c(k) is a constant that depends only on k. The initial condition (t(x, — oo,k) = (|35p together with the fact 
that the asymptotic behavior of (7[x, (3, k) is uniform in /3, and k € (0, 1], imply that c(k) = 0. This gives the following 
expression for c(/3, n): 

r fs^v.,'^ ,87, 



2 y.oo V 
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C. Large-distance asymptotic behavior of the field- field correlator 



As a reminder, we note that the anyonic field-field correlation function can be written combining Eqs. p3p and (|20p 
as follows: 



(vI/t(xi)*(x2))T - ^B++ix, /?, K)e'^(^-'^''') , 7 = (1 + e^n/^ ■ 

Using in this equation the expressions obtained in the previous Section for the potential -B++ and cr, and going back 
to the original variables related to rescaled ones through x = X12VT/2 and /3 = h/T, we arrive at our main result for 
the large-distance asymptotic behavior of the field-field correlator 



i;i2\/TA+ 



0{e 



-2SA+\/T2;l2^ 



(88) 



where 



Ck 



,e"'^^/^a-2(A+) 



(89) 



and are the zeros (|6ip in the upper half-plane. In the asymptotics the functions C(/3, k), c{(3,k), and a{X), 
are defined, respectively, by Eqs. ([4]) or ((87)) . and ((53)) . and the limits of summation over A; are given by ((72)l and 
(|73|1 . The leading part of the asymptotics is the term A: — 0, but as wc approach the free fcrmionic point k 1, the 
A; = — 1 term also becomes relevant. 

VI. ANALYSIS OF THE RESULTS 

As the last step, we check the validity of our main result in the appropriate limits. When the statistics parameter 
K we should reproduce the results obtained for impenetrable bosons in [1, [^, [s^ . Even though in the large- 
distance analysis performed above we have not made any distinction between the cases of negative and positive 
chemical potential (or, equivalently, (3), in what follows, we will see that in the bosonic limit, the asy mptotic behavior 
of the corrrclators is fundamentally different in the two regions, as one would expect from (sl. [sl. [33|. In the fermionic 
limit K 1, our result for the field-field correlator of the anyons should reduce to the correlator of the free fermions. 
At low temperatures, f3 — *■ 00, and positive chemical potential, the system becomes critical, and our result should 
agree with the predictions of the conformal field theory [l^, . 

A. The bosonic limit 

In the bosonic limit k — > 0, Eq. ^ gives in the case of positive and negative chemical potential, respectively, 

e^"-" + 1 



1 f+°° 
lim C(/3, k) = - In 
TT ./ „ 



3^-/3 - 1 



dX + 2iy/j3, /3>0, 



(90) 



and 



1 r+°° fe^^-f^ + A 

limC(/3,.)^-/ In -3;,^^ dA, /3 < . 



(91) 



Also, we have from Eq. ((62 



Introducing the function C(/3): 



A+ = V^, /3>0; A+ = 



/3 < 0. 



(92) 



1 r+°° 

C(/?) = - / In 

TT 



+ 1 



dA, P^h/T, 



(93) 
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and using Eq. (|88p reduced according to the previous equations of this Section, we obtain the following asymptotic 
behavior of the correlation function at negative chemical potential 



(*t(a;i)*(x-2))5 



T 



2^\h\ 



exp 



a{h/T) 



h/T 



-x^2[^C{h/T) + ^\\ 



(94) 



/i < , a;i2 = (xi - X2) ^ oo , 

where 

and we have used (|87p for c(/3, k). Similarly, for positive chemical potential we have 

(*t(a;i)^(x2))T=^e— /l>0, Xl2EE(.Tl-X2)-> 00. (95) 

Both asymptotics ([M|) and ([55]) agree with the result for impenetrable bosons obtained in [3, |4l| . In the case of 
positive chemical potential we made this comparison for the exponential terms only. When /t = and the chemical 
potential (3 is positive, the function go (A) (|54p possesses an index, which means that we cannot transform the RHP 
as in Section flVBI Also, as can be seen from Appendix lEl dC{P,K)/dl3 which appears in ([57)1 becomes divergent 
for K = and /3 ^ 0. This means that the rigorous calculation of the pre-exponential factor for positive chemical 
potential and k ^ would require a more sophisticated approach which we did not attempt in this work. 



B. The fermionic limit 



In the model ([8]) with c oo we use to describe the impenetrable anyons, in the limit of fermionic statistics 
parameter /t — *■ 1, the anyonic field- field correlator should coincide with that of free fermions. To see that Eq. ([55]) 
does indeed reduce to the free-fermionic expression, we note that Eqs. ([1]) - (01) show that for k = 1, 

C(/3,k) = 0, c(/3,k)=0. 

This means that, in this case, a{X) ~ 1, and the leading terms in the correlator ((88)) take the form: 

{^H^i)^{^2))t^^ J2 . (96) 

fc=-l,0 

Taking into account that for k = 1, as follows from Eqs. ([5^ and (|63p . A^j" = — (A^j^)* = (/3 + i7r)^/^, one can see that 
Eq. ([96]) coincides exactly with the asymptotics of the field correlator for free fermions, and reduces to 

i^H^.M^.))T^'^er-^^Y.^^^ (97) 



where 



1/2 _ / \ 1/2 



a= (y/h^ + n^T^-h) /V2, b= ( \/ + n^T^ + h) /V2. (98) 



C. The conformal limit 



For positive chemical potential and low temperatures, /3 ^ oo, the system is conformal. The behavior of C(/3, k) 
in this limit is studied in Appendix [El while the zeros ([S^ and reduce to 
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To make the connection with the conformal field-theory results obtained in [lO|, we use the original variables 
X = {xi — X2)VT/'2,, f3 = h/T and the fact that in the notations used in this work, the Fermi momentum and velocity 
are, respectively, kp ~ ^/h and vp ~ Then, using Eqs. ([55]) and ([M)) . we obtain the following result for the 

asymptotics for the field-field correlator at low temperatures. The accuracy of our calculation of the asymptotics (|88p . 
indicated by the last term in this equation, and leading to the conditions ([7^ and ([75)) . implies that for < k < 2/3 
we can keep only one term in the asymptotic expansion of the correlator: 

For 2/3 < K < 1, however, the two terms in the expansion are legitimate, giving 

{¥ixi)^{x2))T ^ coe"""^('^+^)e"i^'=^'' + c_ie""^'^ H^'^T + i]^^x^2kH^-2) _ (^g^) 
The conformal result obtained in ([Tol. [T§|) is 

{^Hx,)^ix2))T^ J2 i?(g)e-^-^['^''+'^" + ^+'(''+^)1e--'=-(2rf+«). (102) 

Q={N±,d} 

One can see that the leading terms in the sum (|102p . which correspond to Q = 0, 0, and Q = 0, 0, — 1 , are identical 
(modulo some constants) to (jlOOp and (jlOip . The presence of the second term in Eq. (jlOip (and the term with 
Q = 0,0, —1 in the conformal expansion) explains why close to the fermionic point the correlation function exhibits 
beatings, see also Eq. At T = 0, this feature of the correlation function was noticed and explained by Calabrese 

and Mintchev [ll|. 



VII. CONCLUSIONS 



We have computed rigorously the large-distance asymptotic behavior of the field-field correlation functions of 
impenetrable anyons at finite temperature. In the process, we have also obtained a system of differential equations 
which characterize the correlators at any distance. Our result agrees with the predictions of the conformal field theory 
at low temperatures, and describes a transition in the behavior of the correlators between the impenetrable bosons 
realized for the statistical parameter k = 0, and free fermions at k = 1. Variation of the asymptotic behavior of 
the field-field correlation function with the statistics parameter k illustrates the role of the exchange statistics in 
one dimension. This variation contradicts the intuitive notion that exchange statistics is irrelevant in ID systems of 
impenetrable particles, because they can not be exchanged in ID geometry. 

The next important step in the analysis of the ID anyons would be to extend the results of this work to the de- 
scription of the time dependence of the field-field correlators. The resulting time-, space-, and temperature-dependent 
correlation functions could be calculated along the lines of Ref. [3l|, which treated the case of impenetrable bosons. 
We will attempt to do this in a future publication. 



APPENDIX A: SHORT-DISTANCE ASYMPTOTICS 



In this appendix, we obtain the short-distance asymptotics of the potentials and and as a byproduct, 

compute the short-distance asymptotics of the field correlator. First, we need to express the potentials in the form 

more amenable for short-distance computations. We start with Multiplying both sides of the integral equation 

([T7)) which defines /+(A) by 7e_(A) one gets 

+00 Y _ ^-2i{X-fj.)x 



7i?(A)-f7^(A) / ^TTTT — 7/+(M)0%)e-^^"rfM■ 
2^(A-pl) 



This equation shows that B^ can be written as 



B+^{x,p,j) = I siX)dX, (Al) 
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where the function s(A) solves the following integral equation 



s{X) - I ^TTTT ^s{^^)d^l = ■ (A2) 



In a similar fashion, we get 

"+00 



i?++(.T,/3,7)= / e'^^-siX)dX, (A3) 



—00 



where s(A) is the solution of the same integral equation (|A2[) . 
For small x, the solution of (jA2[) can be expanded as 



s(A) = s(A,/3,7)=E«'=(^'/5,7)x', 



k=0 

where Sk are defined by the following recursion relations 

soiX) = 7^(A) , s„(A) = ,so(A) ^ \ ^ _ , , . / (m - A)™-'^-i5,.(/.)dM , m > 1 . (A4) 

k=o 

Defining 



Pi{fi,^)^^l X'^{X)dX, 

'J —00 

i.e., = for odd and using (|A4p in (|Aip and (|A3p . we obtain the short-distance asymptotics for the potentials 
B++{x, (3, 7) = /3o + P^o^ + {Pi - 2/32) x2 + (^/34 _ !/3o/32^ + 0{x^) , 
B+_(.T, (3, 7) = /3o + /^o'a; + /3o'a;' + (^P^ - 1^2^ + 0{x^) . 

These relations give us the short-distance asymptotics of the field correlator. Combining them with Eqs. ([33]) and 
(133 







5(x,/3,7) = S++(x,/3,7)e"^"'^'^^|^=(i+e-»)/., a(x,/?,7) = " / B+.{y, (3,-i)dy , 

we find first 

a{x, (3, 7) = -Pox - ^-Plx^ - ^Plx"" + 0{x^) , 



and then 



/3, 7) = /3o ( I - 2|^x2 + ^/32x3 ) + 0{x^) . 



In the original variables x — {xi — X2)\/T/2 > 0, /3 = /i/r, A ^ X/VT, this result gives for the correlator (fT3)) 

(Vl't(xi)vl/(X2)>T = ^ (1 - ^(^1 - ^2)' + 7^(^i - ^2)') + O ((xi - X2)^) , 7 = (1 + e'^n/^ , (A5) 



where 



1 /•+°° dA ^ 1 f+°° X^dX 



arc the particle and the kinetic energy density, respectively. 



21 



APPENDIX B: LOW DENSITY EXPANSIONS 

As usual, the low-density limit is reached when the chemical potential is such that /3 = h/T — > — oo. In our rescaled 
variables, the density of impenetrable anyons is given by 

^ y/T f+°° dX 



27r 7-00 1 + ' 

so that _D ^ for f3 —oo. In what follows, it will be convenient to use the variable 

C = -e'', C^O for /3->-oo. 
In order to obtain the low-density expansions for the potentials 

oo oo 
A;=l k=l 

we again use Eqs. (jAl[) and (|A3p . and the integral equation (|A2p which has the form suitable for interation expansion 
in density. In terms of the Fermi weight can be represented as 



k=l 



Expanding also s(A): 

oo 

s(A) =^C"sfc(A,a;), 

fc=i 

we obtain from Eq. (jA2p the following recursion relations for the "coefficients" Sk- 
si(A) = -76-^', 

Sfc(A,x) = e Sfc_i(A, x) - 76 / — — — — Sk^i{fi, x)dfi , k>2. 

The first terms of the expansions of the potentials obtained from these recursion relations and Eqs. (|Aip and (|A3p 
are: 

B+_(x,C,«;) = -7V^C+ (-7^ + 7'^^'e--'d.Ti^C' + 0(C3), 

B++{x, C, = -7x/^e--'C + (-^yfe--^ + 7'7re--' e-^^'^+'^^^dx,^ ^ + O(C^) . (Bl) 

Similarly to the short-distance expansions, Eq. (|B1[) gives 

a(a;, /3, 7) = -7V^.Te'^ + ©(e^'^) , 
and the correlator in the rescaled variables: 



In the original variables, this result for the correlator is: 

{¥{xi)^{x2))t = De-^^^i-^^)'/^ (B2) 
It is valid as long as we can neglect the Ole^^) terms, i.e. for T[xi — .12)^ <C \h\/T . 
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APPENDIX C: SOLVABILITY OF THE MATRIX RIEMANN-HILBERT PROBLEM 

As we have shown in Section IIV A[ the matrix RH problem ((42|) is equivalent to the system of nonsingular integral 
equations for functions /±(A): 

/ + 00 
i^T(A,/i)/±(/i)dAi-e±(A), 
-oo 

with the kernel (jl5[) . This means that the RH problem has a unique solution whenever this system of integral equations 
of Frcdholm type has a unique solution. To analyze these equations, we fix [3 and k, leaving the kernel Kt{X, fJ-) a 
function of coordinate x. Let D be an open connected subset of the complex plane and C{Ti.) - the space of operators 
acting on a separable Hilbert space Ti.. Consider the function 

fix) -.D^CiH), 

which for each x in D gives the integral operator with kernel AV(A, fi). Then for each x in the finite strip 

< a <^x <b, 5x<e, 
f{x) is analytic operator-valued fimction. The kernel Kt{X,^) also satisfies the estimate 

/+00 /•+00 
/ |AT(A,/i)|2dA d/^ < C&^ 
'OO J — oo 

where C is a constant, which means that f{x) is compact for each x & D (sec Thm. VI. 23 of [s^). Under these 
conditions, we can apply the analytic Fredholm theorem. 

Theorem C.l (Thm VI. 14 of [soj). Let D be an open connected subset of C Let f : D ^ '^(TY) be an analytic 
operator-valued function such that for each z G D, f[z) is compact. Then, either 

a) (/ — f(z))^^ exists for no z E D; or 

b) (/ — f(z))^^ exists for all z G D\S , where S is a discrete subset of D (i.e. a set which has no limit points in 
D). Ln this case (I — f{z))^^ is meromorphic in D, analytic in D\S, the residues at the poles are finite rank 
operators, and if z G S , then equation f{z)tl) = has a nonzero solution in Ti.. 

As a consequence of this theorem, we have to prove that for at least one point in the strip D the integral equations 
have a unique solution. But this is definitely true for small x, where the Liouville-Neumann scries is convergent. 
Thus, we have shown that the matrix RH problem has a unique solution except for a countable set of values of a;„, 
which we will denote hy X = {a;„}. 



APPENDIX D: SCALAR RIEMANN-HILBERT PROBLEM 



This Appendix provides the basic information on the scalar Riemann-Hilbert problem used in Sec. IIV Bl The general 
problem of this type for the semi-plane is formulated as follows [i^l- Consider two functions g(A) and r(A) defined on 
the real axis, with g{\) nonvanishing. Both are assumed to satisfy the Holder condition: \g{\i)— g{\2)\ < C|Ai — A2|'^, 
and similarly for r(A), with some power k: < fc < 1. One needs to find the function a(A), or a(A), which is analytic 
separately in the upper and lower half-plane, with the boundary values on the real axis satisfying the conditions: 

q:_(A) = a+(A)g(A) , A e M homogeneous problem, (Dl) 

or 

a^{X) ^ a+{X)g{\) + r{X) , AgM inhomogeneous problem. (D2) 

For the purposes of this work, we will assume also the normalization condition a(oo) = a{oo) ~ 1. The considerations 
presented below can be extended to the more general case of a simply-connected closed contour in the complex plane 
- see HO]. 



23 



1. The homogeneous problem 

We need to distinguish three cases depending on the index x(ff) = (l/27r)Var[_oo ^ooiarg g{X) of the function g{X). 
If X = 0, the RH problem with the normahzation condition is uniquely solvable. If the index is positive, x > 0, the 
problem has x + 1 linearly independent solutions, whereas the problem has no solution for x < 0. For x = 0, which 
is the situation most important for the present discussion, the solution of the RH problem (jDl[) is given by 



a{X)=exp\~^ r i^^^rf/iS>, AeC/R. (D3) 



27ri y_oo M - A 

As in the matrix case, it is straightforward to show that the scalar RH problem (jPip is equivalent to the singular 
integral equation 

a+{X) = l + — — dfi, AeK. 

Ztti J_oo jj, — X — tu 

2. The inhomogeneous problem 

Similarly to the homogeneous problem, the solution of the inhomogeneous RH problem (jD2|l with the normalization 
condition is unique for x = 0, which is the situation of interest for the present discussion. The solution can be obtained 
from the solution of the homogeneous problem with the same g{X). If a(A) solves (IDip . then g{X) = a;_(A)/Q;+(A) 
for A e K, and (|D2p can be written as 

a+{X) <5_(A) _ r(A) 



a+(A) a_(A) a_(A) ' 

The functions Q;+(A)/a+(A) and a_(A)/a__(A) are the boundary values of the function a(A)/a(A) which is analytic 
in the complex plane minus the real axis and approaches 1 at infinity due to imposed normalization conditions 
a{oo) = a{oo) = 1. Using the properties of the Cauchy integral, one obtains from this: 



a(A) 27riJ_^ a_(/i)(^-A) 

Thus, solution of the inhomogeneous scalar RH problem (|D2|) is given by 



5(A).a(A)(l-^£"-^*l-^.,), AcC/«^ (D4) 

where a (A) is the solution of the homogeneous problem (|Dip . The singular integral equation equivalent to the 
inhomogeneous RH problem is 

a+(A=l + — / — dn-— -d^l, AeR. (D5 

2tti J_^ fi — X — lO 2^1 J_^ II — X — lO 

APPENDIX E: ANALYSIS OF C(/3, k) 

In this appendix, we study the behavior of the function C{f3,K), which enters the asymptotics (f88|) of the field 
correlator, and is defined by Eq. Q with condition for large and small (3. We start with /3 < 0. In this case, the 
expansion of the logarithms in Eq. ([3]) according to ln(l — z) ~ — Sn=i ^"/"- foi' l-^l < 1; gives: 

C{(3, n) ^ - N ^ + dX , 

^ J-^ .^1 V " n ) 

\ / ^_-[^'jn+lg-n|;3| gm7rKg-n|/3| 



n=l 
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Therefore, the leading terms for large and small |/?| are, respectively, 

'\i3\ p-m 



C(/3, k) = — ^(1 + costtk) + i — 



■ sinTTK , p —> — oo , 



and 



E 



-1)"+^ + COSnTTKN 1 
+ 1—F= 



l3/2 



E 



/TT ' 71' 

n— 1 



3/2 



0. 



(El) 



(E2) 



When /3 is small, its sign is irrelevant and, as we see explicitly below, the last equation holds both for negative and 
positive (3. 

For /3 > 0, we can transform the logarithms in Eq. ([3|) so that the same expansion is applicable, and get 



V0 

2 

TT 



n=l 



dX 



(E3) 



When (3 is large, the formulae 
Jo 

simplify this expression into 



O 



2nV^ / ' 

(—1)"+-^ costittkA 



dX 



O 



^ ' n— 1 



This expression can be transformed finally using the formulae (0.234) and (1.443) of [43, SfeLi (^1)"^^/^^^ = 7'''^/12 



and cosriTTK/n^ = tt'^ B2{k/2), where 



1/6 is the second Bernoulli polynomial. This gives 



C(/3,«) = -^(^y-A^+0~2zV^(K-l), (/3^oo). 
Equation (|E3p also shows that for /? ^ +0, C(/3, k) is given by the same Eq. (jE2p as for /3 ^ —0. 
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